Independent variables of pH: Ten Knights of the Hydrogen Ion Kingdom-Part I. A prospective observational study

CO2, HCO3, SID, and total weak acids have been defined as pH’s independent variables. However, according to Gamble, HCO3 should be equal to the difference between the sum of cations and the sum of anions besides HCO3. Therefore, if this mathematical expression is substituted for HCO3 in the Henderson-Hasselbalch equation, all independent variables of pH can be demonstrated. Our aim is to test this theory in this study. This prospective observational study was conducted between 2019 and 2020. All admitted patients to the intensive care unit who were >18 years old were included. Demographic data, blood gas parameters, albumin, magnesium, and inorganic phosphorus levels, and outcomes were recorded twice (at admission and at the 24th hour). The multivariate linear regression model was used to determine pH’s independent variables. In the multivariate linear regression model, pH was significantly increased by each unit increase in Na, K, Ca, and Mg (mmol L-1). In contrast, pH was significantly decreased by each unit increase in CO2, Cl, lactate, albumin (g dL-1), inorganic phosphorus (mg dL-1), and the strong ion gap. Ten independent variables can accurately predict the changes in pH. For this reason, all ten independent variables should be separately evaluated when interpreting the acid-base status. With this understanding, all algorithms regarding acid-base evaluation may become unnecessary.


Introduction
Independent variables of blood pH must be far more than defined independent variables in time.The essential truth is that hydrogen ion concentration ([H]) is the only determiner of acidity or alkalinization in a solution [1].Over the years, only HCO 3 and CO 2 were used for the acid-base evaluation.From the 1940s to 1970s, many tried to solve the metabolic acid-base abnormalities on the relationships between HCO 3 , CO 2 , the anion gap (AG), the buffer base concept, and standard base excess (SBE) [1][2][3][4][5].However, the most crucial theory in these years was Gamble's law of electroneutrality: "The sum of cations equals the sum of anions in the plasma" [6].In the 1980s, Stewart proposed three independent variables for pH by using the law of electroneutrality: strong ion difference (SID), total weak acids (A tot ), and CO 2 [7].This was the first theory that mentioned independent variables of [H + ].In the same decade, corrected AG (AG corr ) was added to the traditional approach algorithm, and then Fencl derived his algorithm based on the Stewart approach [8,9].This millennium, O'Dell and Story used partitioning base-excess models to solve metabolic acid-base disorders by separating the effects of Na-Cl difference, albumin, lactate, and unmeasured anions on the SBE [10,11].However, none of these approaches have mentioned all the independent variables of [H + ].Yet, according to the Gamble equation, HCO 3 should be equal to the difference between all cations' ionic charges (Na, K, Ca, Mg) and all anions' ionic charges other than HCO 3 (Cl, albumin, inorganic phosphorus (P i ) , and unmeasured anions (UA)).In 2005, Schu ¨ck finally derived the new Henderson-Hasselbalch equation by substituting Stewart's independent variables for HCO 3 but there was no mention of the independent variables of [H + ] in his study [12].Unexpectedly, he used Stewart's independent variables and AG corr together to define the metabolic components.In this study, we sought to determine all independent variables of [H + ] using the Gamble and Henderson-Hasselbalch equations and we hypothesized that there were ten independent variables of pH.

Formulas
1. Mg and P i values were converted from milligrams to millimoles by using the formulas below [16,17]: Mg (mmol L -1 ) = Mg (mg dL -1 ) x 0.41152 P i (mmol L -1 ) = P i (mg L -1 ) x 0.323 2. [Alb -] and [P i -], A TOT , SIG, AG, and AG corr were calculated by using the formulas below [12,18]: Blood samples were taken from the invasive arterial cannula of all patients at two different times: at the ICU admission and the 24 th hour of the ICU period.For blood gas samples, syringes with dry heparin were used, and the samples were tested by a blood gas machine in the ICU.As for biochemical parameters, syringes without heparin and tubes were used, and these samples were tested in the biochemistry laboratory.Blood gas data were acquired using an ABL 800 (Radiometer, Denmark, Copenhagen) blood gas device, which employed ion-selective electrodes.Alb, Mg, and P i were acquired using a Cobas C 303 device (Roche, Rotkreuz, Switzerland).

Theory
1.According to the Henderson-Hasselbalch equation, there are two independent variables of pH: 2. Gamble JL claimed that the sum of cations equaled the sum of anions in the plasma: [6] Na 3. According to Gamble's equation, HCO 3 can be written as below.4.This mathematical expression can be written as a form of the Henderson-Hasselbalch equation instead of in terms of HCO 3 .Since UA's ionic charge is calculated using SIG [14], the Henderson-Hasselbalch equation can be revised as below: We tested whether all parameters in the last Henderson-Hasselbalch equation were independent variables of pH.

Statistical analysis
Descriptive data are presented as mean±sd, median (quartiles), and percentages.The Kolmogorov-Smirnov test was used to detect normality.Multivariate linear regression models were used to detect independent variables of pH by adding P a CO 2 , HCO 3 , Na, Cl, K, Ca, Mg (mmol L -1 ), lactate, Alb (g dL -1 ), P i (mg dL -1 ), and SIG.The Enter method was used to build the model.Furthermore, two multivariate linear regression models were used for the traditional and Stewart approaches.To determine the best model for predicting the pH, leave-one-out cross-validation (LOOCV) scores were used.The correlations and measures of agreements between AG corr and SIG were analyzed with Pearson correlation and the Kappa test, respectively.The estimated power of this study was calculated as 0.99 using the F test for the regression model (effect size f 2 = 0.35, α = 0.05, total sample size = 410, and the number of predictors = 10, GPower 3.1.9.4).SPSS version 29 was used for all statistical analyses.p<0.05 was accepted as significant.

The traditional approach
This approach is based on the H 2 CO 3 -HCO 3 buffer system [1,19,20].CO 2 and HCO 3 are deemed independent variables.Based on this, acidosis is defined as an increase in CO 2 or a decrease in HCO 3 , and vice versa for alkalosis [19,20].However, CO 2 is a measured variable, whereas HCO 3 is a calculated variable.Moreover, since HCO 3 alone is not enough to detect metabolic components, other calculated variables, such as SBE and AG, must be added to the evaluation [21].These calculated variables also include other variables in their formulas, including HCO 3 [21].In addition, AG should be corrected for albumin at least [22].To detect whether there are compensation disorders and mixed disorders, rule-of-thumb equations and expected (standard) HCO 3 are calculated again and assessed; then, a decision about acid-base status is made [19,21,23].For these reasons, the traditional approach has substantial limitations.Firstly, HCO 3 is a calculated variable and cannot be an independent variable for pH.Therefore, other parameters derived from HCO 3 , which are SBE and AG, cannot be independent variables as well.Furthermore, AG should be fully corrected, but if it is fully corrected, it becomes exactly the SIG of Stewart's approach.For this reason, the measures of agreement between albumincorrected AG and SIG may be weak (Table 3).Therefore, metabolic acid-base disturbances are tried to be understood by using these dependent variables.Secondly, it is perceived as if the H 2 CO 3 -HCO 3 buffer system is the only system for [H] production.Thus, other probable metabolic variables that affect water dissociation are ignored or overlooked.Further, according to this approach, respiratory disturbance (CO 2 ) is compensated for by only metabolic components (HCO 3 ) and vice versa [19,23].For this reason, compensations or different effects in the metabolic components remain unspecified (e.g., hypochloremic alkalosis and hyperlactatemia under normal CO 2 levels).Thirdly, algorithms created based on these facts are still complicated and confusing.According to many textbooks, it would seem like there are no acid-base disturbances if pH, CO 2 , and HCO 3 are normal [15,24].We believe that all these reasons may explain why the traditional approach to predicting pH is the weakest (Table 2 and Fig 3).

Stewart approach and then: New calculated independent variables and corrections
Stewart's approach creates a solution to metabolic acid-base disturbances by using the electroneutrality law, which Gamble summarizes [6,7,18].In this approach, the independent variables of pH are CO 2 , SID, and A TOT .SIG was also added to Stewart's independent variables [14,18].SID, A TOT , and SIG are calculated variables, whereas CO 2 is a measured variable in this approach.In this situation, SID and A TOT also relate to their components in both equations [18].These relationships can cause some problems.For instance, SID can be normal or high in hyperlactatemia, and A TOT can also be normal or low in hyperphosphatemia.To explain these unexpected situations, we must look deeper into Stewart's equations.Yet, by its definition, an independent variable should have a clear and indisputable effect and thus should not cause these unexpected situations.According to our results, all SID and A TOT components are independent variables of pH (Table 2).Therefore, creating new calculated variables by combining these independent variables is unnecessary.In these circumstances, the only exception is SIG.Interestingly, we found that SIG was an independent variable in this study, although it was a calculated variable (Table 2) [14,18].Additionally, it was the fourth most important independent variable of pH (Fig 2).Despite these limitations, the prediction of pH using Stewart's approach is better than using the traditional approach (Table 2 and Fig 3).In this millennium, the partitioned BE model and chloride corrections have been devised, referring to solving metabolic disturbances [10,25].Partitioned BE is a useful method at the bedside, but the relationship of pH with its components has yet to be mentioned.On the other hand, chloride corrections, put forth by Fencl, correct the observed serum Cl level with the serum observed Na level [9].If the newly calculated Cl level is not between the defined normal limits, it is accepted as hypo-or hyperchloremia.This point of view does not make sense because the observed serum Cl level is measured and is the actual value for serum Cl, whereas corrections are calculated and imaginary values for serum Cl.Measured Cl is a real molecule, and this real molecule has a real ionic charge.

What is the meaning of the existence of ten independent variables?
They have always been there; however, we prefer not to see them as individuals.
This perspective, which is the best at predicting pH (

Water resolvents (electrolytes and lactate)
According to Stewart's approach, the difference between strong ions affects water dissociation.Combining this knowledge with our results, we can conclude that deviations from normal serum values of each electrolyte and lactate affect water dissociation only (Fig 4).Among these, Na, Cl, and lactate ionic charges have substantial effects on [H] separately (Fig 2 ).

Buffers (albumin and ionic phosphorus)
Albumin and P i are named as weak acids in Stewart's approach [7,18].They do not affect [H] on water dissociation but bind hydrogen ions when hydrogen ions increase in the plasma and vice versa (Fig 4).Indeed, these buffers' ionic charges are affected by the pH in accordance with their ionic charge formula [12,18], but at the same time, they are also independent variables of pH (Table 2).This indicates that even though albumin and Pi serum levels do not change, their ionic charges vary when pH increases or decreases.Another meaning is that albumin and P i continuously correct the changes in pH due to the effects of other independent variables.Albumin may do this by changing its structure [26].On the other hand, decreases and increases in albumin result in alkalosis and acidosis, respectively [9].This situation can only be explained by the existence of an inverse relationship between albumin level and its affinity to hydrogen ions.Albumin's affinity to hydrogen ions increases in hypoalbuminemia and decreases in hypoalbuminemia.Interestingly, this relationship is the exact opposite of that between hydrogen ions and hemoglobin, another essential plasma buffer.It seems that albumin works to protect the oxygen-binding capacity of hemoglobin because the most crucial duty of hemoglobin is oxygen transportation.

Hydrogen ion producers (CO 2 : An indisputable independent variable of all approaches)
CO 2 is neither a water resolvent nor a buffer in the plasma.It is only a hydrogen producer in the H 2 CO 3 reaction (Fig 4).For this reason, it is an independent variable of pH and has the most important effect on it (Fig 2).However, this effect should involve the Hb buffer effect because CO 2 transporting occurs via the same reaction in erythrocytes.On the other side, although HCO 3 seems to be a buffer in this reaction, its buffer effect was not enough to make it an independent variable of pH in our results (when we added HCO 3 to Model III, the pvalue for HCO 3 was 0.125).HCO 3 and hydrogen ions are already end-products of this reaction and are associated with CO 2 .Therefore, the outlook on this reaction should be changed.
Unmeasured anions (beta-hydroxybutyrate, D-lactate) and cations (Fe, Li, globulins) Unmeasured anions/cations are undoubtedly important independent variables of pH (Table 2 and Fig 2).However, we still have to calculate their effect using the SIG formula, which may be a limitation [14,18].We believe that our new approach in Part II can logically solve this limitation.UA may have different mechanisms to affect [H].Therefore, we have avoided categorizing their effect mechanisms (Fig 4).The study's most important limitation is that it is a single-center study.On the other hand, as new metabolic parameters are started to be measured routinely over time, these independent parameters may increase.Thus, this topic may be verified and improved by multicentric prospective studies in the future.

Conclusions
pH has ten independent variables, and we reach the best prediction for pH when using all of them.Therefore, each of them should be separately evaluated when interpreting acid-base status.Each independent variable has an individual acid-base effect on top of its clinical importance.For this reason, it is unnecessary to create calculated parameters from them.On the basis of our results, we argue that blood gas machines have to measure each independent variable because it is their only purpose of existence.Evaluating these independent variables separately will yield more authentic results, rendering all complicated algorithms unnecessary.

Table 2 and
Fig 3), directly links each independent variable to[H].At this point, it is essential to discuss how and in which direction they affect [H].